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Abstract T violation has previously been shown to induce destructive interference between different 
paths that the universe can take through time and leads to a new quantum equation of motion called 
bievolution. Here we examine further details of the interference and clarify the conditions needed for 
the bievolution equation. 
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1 Introduction 

Parity inversion (P), charge conjugation (C) and time reversal (T) are fundamental operations in 
quantum physics. The discoveries that Nature is not invariant under P, CP and T stand as pillars of 
20th century physics that promise far reaching consequences. The implications of the violation of CP 
invariance for quantum decoherence, Bell inequalities, direction of time, particle decay, uncertainty 
relations and entropy generation has been explored over the last few decades [IH25]. Less attention has 
been given to the violation of T invariance, which is implied by the violation of CP invariance under 
the CPT theorem and is also observed directly [3S1I3Z]- Nevertheless, this fundamental time asymmetry 
has been linked to the direction of time [38l[39]. In addition, I have recently shown that the violation 
of T invariance in neutral meson systems has the potential to affect the nature of time on a large 
scale 130] far beyond that previously imagined [33]. 

The large scale effects of the violation of T invariance (hereafter referred to as simply T violation) 
are seen to arise when modelling the universe as a closed system with no predehned direction of 
time [40] . T violation implies that there are actually two versions of the Hamiltonian, Hp and Hp, one 
associated with evolution in each direction of time. To avoid prejudicing one direction of time over the 
other, the time evolution of the system is modelled as a quantum walk comprising a superposition of all 
possible paths that zigzag forwards and backwards in time. The size of each step in time is fixed at the 
Planck time. Interference due to T violation is found to eliminate most of the possible paths leaving a 
superposition of just two paths: one comprises continuous evolution under Hp in the forwards directions 
of time and the other comprises continuous evolution under Hp in the backwards directions of time. 
This result has been called the bievolution due to the evolution under two Hamiltonians. Standard 
quantum mechanics is recovered in each path. However subtle features of the interference that were 
not previously appreciated have come to light. The purpose of the present paper is to reveal those 
subtle features and explore their repercussions. We shall see that the regime under which bievolution 
is valid is more restricted than previously thought. However, we shall also find that by reducing the 
size of the time step appropriately, the newly-found restrictions are avoided. 
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We begin with a brief review in Sect. 2 and then draw out the subtle details of the interference in 
Sect. 3. In Sect. 4 we show how they impose a limitation on the range of validity of the bievolution 
equation and then we examine a resolution of the problem in Sect. 5. We end with a discussion in Sect. 

6 . 


2 Review of earlier work 


We briefly review the results presented in Ref. m where full details can be found. Our analysis 
applies to a model universe that is composed of matter and fields on a non-relativistic background 
of space and time. This model is imagined to be comparable in scale with the visible portion of the 
actual physical universe. One could subdivide the model universe into parts which act as spatial and 
temporal references for the remainder. We assume, however, that the model universe itself is not a 
subdivision of something larger, and so there is no physical system external to it. Thus, there is no 
external object that could be used as a reference for an origin or direction of space for the universe 
and, likewise, there is no external clock or memory device that could be used as a reference for time. 
In particular, this means we have no basis for favouring one direction of time evolution of the universe 
over the other. However, maintaining a description that is unbiased with respect to the direction can 
become cumbersome and so, for convenience, we arbitrarily label the time evolution in the directions 
of the positive and negative time axes as “forward” and “backward”, respectively. 

Let the evolution of the universe in the forward direction over the time interval r be given as 
IV'J’(r)) = Uf{t)\iIjo) where 

Uf{t) = exp(—. (I) 

Here and elsewhere we set h = 1 for convenience. Similarly, let the evolution in the backward direction 
by the same time interval be |'0 _b(t)) = IJb{t)\'4’o) where 

Ub{t) = exp{iTHB) ( 2 ) 

with Hb = THfT~^, Ub{t) = TUf{t)T~^ and TiT~^ = —i where T is Wigner’s time reversal 
operator [45) . In other words, Hf and Hb are the generators of time translations in the positive—t 
and negative—t directions, respectively. 

If observers within the universe were to consistently find experimental evidence of one version of the 
Hamiltonian, say JLf, it would imply that the time evolution is in the corresponding direction, i.e. in 
the positive-t direction. But as our analysis needs to be unbiased with respect to the direction of time, 
it must include both versions of the Hamiltonian on an equal footing. An analogous situation occurs 
in the double slit experiment: if there is no information regarding which slit the particle goes through, 
we must add the probability amplitudes for the particle going through each slit. With this in mind, 
consider the expressions (^|C/f(t)|^o) and {4>\UB{T)\'ij}Q) which represent the probability amplitudes 
for the universe to evolve from the state |'0o) to the state \(j)) via two different paths, each of which 
corresponds to a different direction of time and associated version of the Hamiltonian. As we have no 
reason to favour one direction (or Hamiltonian) over the other, we follow Feynman’s sum over paths 
method and take the total probability amplitude to evolve from I'i/^o) to \(j)) as the sum 

{(I}\Uf{t)\iPo) + {(I)\Ub{t)\iPo) = {(I)\[Uf{t) + C/b(t)]|?/'o) • (3) 

As this result is true for all states |(/)), it represents the universe evolving from |^o) to the state 


mr)) 


Uf{t) + Ub{t) 


IV’o) ■ 


(4) 


We call the process described by Eq. symmetric time evolution over one step in time. Symmetric 
evolution in both directions of time has also been explored by Carroll, Barbour and co-workers [1TH45] 
in different contexts. Here it is simply a result of avoiding any bias in the direction of time. 

The extension to the case where there are many different possible paths is straightforward: the total 
probability amplitude for the universe to evolve from one given state to another is proportional to the 
sum of the probability amplitudes for all possible paths through time between the two states. We must 
be careful, however, to weigh all paths in the sum uniformly. We can do this in a recursive manner by 
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replacing |'(/;o) in Eq. ([3|) with |S'(t)) to give the sum of the amplitudes for evolving from |^'(t)) to \(j)) 
via two different paths in time. The symmetric time evolution of that results is given by 


\^i2T))= UFiT)+UBiT) \WiT)) 


(5) 


which, according to Eq. 0, is the two-step symmetric evolution of \tpo)'. 


m2r)) 


Uf{t) + UsiT) 


2 


IV'o) • 


( 6 ) 


Repeating this argument N times results in the A'^-step symmetric evolution of |^/’o) as 


mNr)) 


IJf(.t) + Ub{t) 


N 


l^o) • 


(7) 


This was shown to be equal to 

N 

\^{NT))=Y,SN-n,n\^^) (8) 

n—0 

where SM-n,n is the sum of (^) terms each of which is a unique ordering oi N — n factors of Ub(j) 
and n factors of Uf{t). 

Any given ordering of factors, such as ■ ■ ■ Uf{t)Uf{t)Ub{t), represents the evolution along a 
path that zigzags in time, going one step backward, two forward, etc. and so the expression {ip\ ■ ■ ■ 
UFi,T)UF{T)UBiT)\'ipo) represents the corresponding probability amplitude for evolving from |^o) to 
\lp) over that path. This means that the expression {^p\SN-n,n\'^Q) represents the total probability 
amplitude for evolving from |^/>o) to \^p) along (^) different paths in time. The evolution along different 
paths can destructively interfere leading to a small, or even zero, total probability amplitude. A clock 
device that is part of the Universe and constructed of matter that obeys T invariance, would evolve in 
time by the same net amount of {N — 2n)r for each path represented by {lp\Sm■ The use of 
T-invariant matter in the construction of the clock device avoids any ambiguity in the length of time 
intervals for different directions of time. 

Reordering Sm,n so that all the Ub{t) factor are to the left of the Uf{t) factors gives 


Sr. 


m si 

UB{mT)UF{nT) E-EE exp 
^=0 A;=0 


{v 1 -\- k^T^\HF. Hb] 


(9) 


to order r^. There are n summations on the right side and [A,i3] is the commutator of A and B. It 
is useful to express the summations in Eq. dH]) in terms of the eigenvalues of the Hermitian operator 
i[HF,HB]- Writing the operator i[HF,HB] in terms of its eigenbasis as 

i[HF.HB] = J \p{X)fl{\)dX , ( 10 ) 

where 11 (A) is an operator with unit trace such that p{X)n (A) is a projection operator that projects onto 
the manifold of eigenstates of i[HFT Hb] with eigenvalue A, and p(A) is a density function representing 
the degeneracy of eigenvalue A with p(X)dX being the number of eigenstates in the interval A, then 


yields [JD] 

Sm.n = UBirrn 

-)Uf(nT) J Irr.,n{r‘^A)p{X)n{X)dX 

(11) 

whertQ 

m 

s 1 



^ ^ ■ 

■ ■ ^ ^ exp [— f (r; H- + (■ + k)z\ 

(12) 


v=0 i=o k=0 


^ Note that /m,n('r^A) here is equivalent to the function 7m,n(A) in Ref. |40| . 
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is an interference function which has been shown to be symmetric with respect to the indices n and 
TO, i.e. Im,n{z) = In,m{z). After some algebraic manipulation [40] we find 


In-t, 


n 

exp[—m(A^ — n)z/2]]^ 
9=1 


sin[(iV + 1 - q)z/2] 
sin(gz/2) 


(13) 


In Ref. [30] the value of r was fixed at the Planck time r « 5 x 10“^'^s and the degeneracy p(A) 
was taken to be Gaussian distributed about a mean of A = 0. It was shown that /Ar_„^„(T^A) was 
sharply peaked at A = 0 for n ^ 0 and n N from which it was argued that the integral in Eq. m 
becomes approximately proportional to the projection operator p(0)7I(0), and so if p(0)7J(0)|^o) = 0 
then the corresponding term SN-n,n is missing from the sum in Eq. An absence of the term 
indicates destructive interference between the associated paths through time. It was concluded that 
for suitable states of the Universe, destructive interference eliminates all paths except for those that 
are either continuously backward (corresponding to n « 0) or continuously forward (corresponding to 
n ~ TV). The result is the bievolution equation 


\mr)) 


Uf{Nt) + Ub{Nt) IV'o) 


(14) 


for sufficiently large values of TV. However, subtle details of the interference function /jv_„,„(r2A) 
that were not previously appreciated place restrictions on the range of validity of this result. In the 
remainder of this paper we explore those details and elucidate their repercussions. 


3 Details of the interference fnnction 


We noted previously that |/Ar_„_„(z)| has a central maximum at z = 0 for which [40] 


lN-n,ni0) = in) ■ (15) 

We shall refer to this as the principle maximum of \lFi-n,n{z)\ and any other maxima as subsidiary 
maxima. We can find other features of the interference function lN-n,n{z) by determining the points 
along the z axis where it takes on particular values. To do this, consider the factors in the iterated 
product of Eq. (fOl) each of which comprises a numerator of the form /num(TV, q, z) = sin[(TV +1 — g)z/2] 
and a denominator of the form fden{q, z) = sin(gz/2). We now deduce three properties as follows. 

(i) The zeroes of lN-n,niz) occur at the points away from the origin where one of the numerators 
is zero, that is, for /num(TV, q,z) = 0 and z ^ 0 which is satisfied by 


2mTT 

N + l-q 


(16) 


for q = 1, 2,..., n and non-zero integer to. 

(ii) The interference function has a modulus of unity at the points where each factor in the iterated 
product has a modulus of unity. This occurs at a nonzero value of z where /num(TV, g, z) = ±fden{q, z) 
for all values of q in the range 1, 2,..., n. Writing this condition as sin[(TV-|- I)z/2 —gz/2] = ± sin(gz/2) 
shows that a set of solutions is given by (TV -|- l)z/2 = tott, or 


2mTT 

N+1 


(17) 


for non-zero integer to. The iterated product in Eq. m can be rearranged to give a product in which 
the numerators are multiplied in reverse order as follows 


n 


n 


/num(-^; Q') 
fdeniq.z) 


n 


n 


/num(TV, 77 . -t- 1 z') 

/den(d; z) 


Each factor on the right side has a modulus of unity for sin[(TV — 77)z/2 -|- gz/2] = ±sin(gz/2), which 
is satisfied by (TV — n)zj2 = tott , i.e. by 

2m'!T 

N — n 


z = 


(18) 
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(iii) The modulus of the interference function becomes maximal at the points where all the numer¬ 
ators satisfy fnum{N,q, z) « ±1 and all the denominators are small, i.e. where fdeniq,z) <C 1 and so 
s\n{qz/2) K. qzjl <C 1 to first order in qzjl. At these points the modulus of the interference function 
is given approximately as 

-‘--‘- 9=1 qz/2 n\ \zJ 

This simple expression will give the approximate magnitude of a maximum in \lM-n,n{z)\ provided we 
can specify the points along the z axis where the maxima occur. These points can be found for n ^ N 
as follows. We restrict our attention to the case where each denominator fdeniq, z) ioi q = 1,2,... ,n 
varies over much less than one cycle with respect to z, i.e. 0 < z <C 47r/n. The condition on the 
numerators implies that sin[(Ai -b 1 — q)z/2] ^ ±1 which, given that n N and so q N, is satisfied 
by 

{N+ l)z/2^ (2m + l)Tr/2 (20) 

for integer m. We shall take first subsidiary maximum along the positive 2 axis as the one that occurs 
after the first zero of Ijv-n,n(z) which, according to Eq. (ITB)) . is for z > 2tt/N. Equation (EUl) then 
gives the position of the m-th subsidiary maximum of \lN-n,n(z)\ along the positive z axis as 


(2m + l)7r 

{N + l) 


( 21 ) 


for positive integer m. Combining this with the condition z Axr/n for the denominators to be small 
implies that Eq. m is valid for the m-th maximum provided 


(2m + l)7r 47r 

(N + l) n 


( 22 ) 


To summarise, this analysis implies that \lN-n,n(z)\ ranges from (i) zero at the points given by 
Eq. (1161) through (ii) unity at the points given by Eq. (dzi) and Eq. dUD and to (iii) subsidiary maxima 
that are bounded by Eq. dig) at the points Eq. (EH- 

The principal maximum is significantly larger than the bound on the subsidiary maxima given by 
Eq. dHI) for values of n that are signihcantly different from 1 and N where N 1. Eigure[T] compares 
the relative heights of the principal and subsidiary maxima for N = 8000 and shows that the subsidiary 
maxima for the cases n = 10 and n = 50 are negligible (and below the resolution of the figure) in 
comparison to the principal maximum. 

Nevertheless the subsidiary maxima play a signihcant role that was not appreciated previously. 
In Fig. [2] we plot the natural logarithm of |/Ar_„^„(z)| for the same cases as in Fig. [TJ Also plotted 
are the bounds (dotted grey curves) given by Eq. (Hi?)) which accurately estimate the magnitudes of 
the subsidiary maxima. Knowing the magnitude of the principle maximum, from Eq. m, and the 
bound on the subsidiary maxima, from Eq. dH, allows us to scale the interference function so that all 
maxima are of equal height using the following scaling function. 


FN,n(z) 


(n) , for a; < ^ 

^ (I)" , otherwise . 


Figure [3] shows the scaled interference function YN-n,n{z) given by 


Y]Sl—n,n(z) 


^N—n,n{z) 
FM ,n(z) 


(23) 


for the same cases as in Fig. [31 Notice that the height of the sixth subsidiary maxima (i.e. for m = 6) 
for the case n = 50 is slightly larger than unity due to the approximate nature of the bound in Eq. m 
used for the scaling function (despite condition Eq. (1221) being satisfied). 

We can estimate the widths of the subsidiary maxima by htting a quadratic function as follows. 
We set 

z = Zm + e , 
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Fig. 1 (Colour online) The modulus of lN-n,n{z) for N = 8000 and n — 1 (red curve), n = 10 (green curve) 
and n = 50 (blue curve). For clarity, each curve has been scaled by dividing |/jv-n,Ti(z)| by {„) to give a value 
of unity at a = 0. The dotted gray line represents the bound on the subsidiary maxima for n — 1 given by 
Eq. (II3. 



Fig. 2 (Colour online) Logarithm of \lN-n,n{z)\ for N = 8000 and n = 1 (red curve), n = 10 (green curve) 
and n = 50 (blue curve). Again, for clarity, each curve has been scaled by dividing ln[|JAr- ra,n( .2)|] by ln[(()^)] 
to give a value of unity at 2 = 0. The dotted gray lines represents the bounds given by Eq. 


where Zm is the position of the m-th maximum given by Eq. (EH), and express the right side of Eq. m 
as a function of e. The quadratic approximation will only be valid over a range of values of e that is 
much less than half the distance between consecutive maxima, i.e. for 


|e| < 


TT 

tvTT ■ 


(24) 


In the limit of large TV, we find that the factors of the iterated product in Eq. (fT51) can be approximated 
for q N and qZm <C 1 by 


sin[(iV + 1 - q){zm + e)/2] ^ ± cos[(A^ + 1 - q)e/2] 
sin{q{zm + e)/2) sin{qZm/2.) 


(25) 


where the sign depends on the value of m and we have made use of Eq. (1241) in the denominator. The 
arguments of the trigonometric functions on the right of Eq. EH are much smaller than tt, and so 
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Fig. 3 (Colour online) The scaled interference function YM-n,n{z) for N = 8000 and n = 1 (red curve), n = IQ 
(green curve) and n = 50 (blue curve). 


expanding them to second order yields 

sin[(fV + 1 - q)[zm + e)/2] 


± 


1 - (TV + 1 - g)2eV8 


sin(g(2;m + e)/2) qZraj^ 

Substituting this result into Eq. (US gives the modulus of the interference function as 

1 - (TV + 1 - g)2e2/8 


(26) 


\^N—n,n{Zm ~ 

9=1 


qZm/‘^ 


which is approximately 


\lN—n,n{,Zm c)| 


n\{zm/‘2.Y 


to second order in e. Evaluating the summation and retaining terms of order {N + 1) or higher then 
gives 


\^N—n,n{,Zm T c)| 



n){N + 1) 
"8 


(27) 


We found using the same technique in Ref. m that the quadratic approximation to the principle 
maximum was given by 


|7v-n.„(e)|^(ir) 


1 - e- 


n{N — n){N + 1) 


24 


(28) 


Figure 0] illustrates the closeness of these approximations for the principle and first subsidiary maxima 
for various values of n. 

The widths of the peaks can be estimated from the points where the right sides of Eqs. (EHl) and 
(EZD are zero. The corresponding values of e are 


^Drin — 


24 


i(iV-n)(7V + l) 


i{N — n){N + 1) 


(29) 


^sub — 


(30) 
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Fig. 4 (Colour online) Quadratic approximations (dotted curves) given by Eqs. (1271) and (12811 to the scaled 
interference function Yjv-n,Ti(z) (solid curves) for the principle and first subsidiary maximum. The plots are for 
N = 8000 and n = 1 (red), n = 10 (green) and n = 50 (blue). Note that different scaling is used for the 2 axis 
for the principal and subsidiary maxima. The subsidiary maximum is centered on z = a where a = 37r/(AC+l) 
and b = 0.00025. 


for the principle and subsidiary maxima, respectively. These values are bounded by 


where k is the lessor of n and N — n. 


^prin ^ 
Csub ^ 


4.9 

VkN 

2.8 

y/kN 


(31) 

(32) 


4 Implications for the bievolution equation 

We have now uncovered sufficient information about the subsidiary maxima to be able to explore their 
implications for the derivation of the bievolution equation Eq. (HI- Combining Eqs. (|H]) and m gives 
an expression for the state \^{Nt)) in terms of the interference function /m,n('r^A) and the density 
function p(A) as 

N 

\mr)) = E Ub[{N — n)T)]Up{nT) 

n—0 

It is important to keep in mind in the following that the state \^{Nt)) is not normalized; in particular, 
to enable a consistent probabilistic interpretation its magnitude needs to be normalized for each value 
of N. In order for Eq. dSSl) to yield the bievolution equation, the terms for values of n that significantly 
different from 0 and N must be negligible in comparison to the remaining terms. The only way for 
this to occur independently of the general nature of the operators Ub{‘)Uf{-) is if the expression in 
square brackets vanishes for these terms. To see how this might happen consider the two functions 
lN-n,n{T^^) = lN-n,n{z) and p(A) = p{zfT^) that appear in the integral of Eq. (1551) . where z = r^A. 
Fig. [S] illustrates the degree to which their product contributes to the integral. The scaled interference 
function Y)v_„_„(z) given by Eq. (1231) is used in the figure to emphasize the role of the subsidiary 
maxima and, as in Ref. m, we assume the eigenvalue density function p(A) is Gaussian and centered 
on A = 0 with a standard deviation of Asd ~ VJ ^ where / represents the fraction of particles 

in the universe that violate T invariance. We also assume that the state j'i/'o) satisfies the nonzero 
eigenvalue condition |40) : 


lN-n,n{'^'^Wdx \tpo) 


(33) 


p(0)77(0)|^o)=0. 


(34) 
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Fig. 5 (Colour online) Sketch comparing the scaled interference function of \YN-n,n{z)\ (red curve) with the 
eigenvalue density function p(zfT‘^) (green curve). 


i.e. it has no projection onto the subspace of eigenstates of i[Hp,HB] with eigenvalue zero. 

First take the case where the position of the m = 1 subsidiary maximum is significantly greater 
than the width of the density function, i.e. where 


Stt 

iV + 1 


> t^Asd • 


(35) 


This places an upper bound on the value of N. In this case, provided 0 ^ n <C fV the integral in 
Eq. (1331) is approximately proportional to the projection operator p(0)77(0) which projects onto the 
subspace with zero eigenvalues. Due to the nonzero eigenvalue condition Eq. (IMl) . the corresponding 
terms in Eq. (1331) are negligible. The bievolution equation then follows, as reported in Ref. |40) . 

Next take the case where the value of N violates the upper bound in Eq. (EH) and, to be specific, 


let 


37r 

N + 1 


t^Asd ■ 


(36) 


Again consider the range 0 n iV and, due to the nonzero eigenvalue condition, note that 
p(O)7T(O)|'0o) = 0. The position of the subsidiary maximum is now centered on z = t^Asd and so 
the expression in square brackets in Eq. dMl) is approximately proportional to p(Xsi))n{XsD)\''Po)^ 
which is nonzero in general. This implies that the corresponding terms in Eq. dMl) are not necessarily 
zero and so the bievolution equation no longer holds. Moreover, further consideration along these lines 
will show that the bievolution equation is not recovered, in general, for larger values of N. 

Thus, Eq. (1351) represents a condition that must be satisfied in order for the bievolution equation 
to hold. Expressing the condition as an upper bound on the total time Nt gives 


Nt <C 


37r 
T Asd 


10-13 s 

~7r 


(37) 


For the right side to be 10 billion years, the fraction / of particles in the visible universe that are 
kaon-like and contribute to T violation needs to be 10“®!. Given it contains an estimated 10®° proton¬ 
like particles, this means that the visible universe needs to contain much less than a mole of kaon-like 
particles in order for the mechanism to explain the direction of time. 


5 Reducing the value of r 

The physical unreasonableness of such a small value of / calls for a review of the approach we have 
employed. Up to now the size of the step in time, r, has been fixed at the Planck time. But Eq. (03 
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can be satisfied for more reasonable values of / if r is allowed to have a much smaller value. Moreover, 
the preceding analysis shows that condition Eq. inevitably fails at some value of N because the 
first subsidiary maximum moves towards the origin of the z axis as the values of N increases. The 
positions of all the subsidiary maxima would remain relatively fixed, however, if the size of r scales 
as l/^N + 1. For example, replacing r in Eq. with c/^/N + 1, where c is a constant, yields the 
constraint on the position of the first subsidiary maximum as 

Stt 

N+1 ^ ’ 

which is satisfied, independently of the value of N, provided c^Asd "C Stt. Correspondingly, the argu¬ 
ment of YM-n,n{z) in Eq. (1^ should be replaced with r^A = c^X/{N + 1) = z/{N + 1). This yields a 
new function of z as follows: 

yN-n,n{z) = yN-n,n{j^) ■ (38) 

Figure |6] verifies that the position of the first subsidiary maximum of \yN-n,n(z)\ converges to 

z = Stt 

for large values of N. As z = c^A, this corresponds to A = . Thus, if yN-n,n{z) is plotted as a 

function of A, the first subsidiary maximum will be fixed at A = Stt/c^. This position can be made far 
beyond the width Asd of the eigenvalue density function p{X) by choosing the value of c accordingly. 
That being the case, the integral in Eq. (I33|) will be approximately proportional to the projection 
operator p(0)J7(0) for 0 <C n ^ iV and, if the nonzero eigenvalue condition Eq. (IM)) holds, the 
corresponding terms in Eq. (13311 will be negligible. Hence, the conditions for the bievolution equation 
Eq. (imi can be satisfied as N increases indefinitely. However, the full analysis of the consequences of 
reducing the value of r in this way is beyond the scope of this work and will be explored elsewhere |46j . 

6 Discussion 

In Ref. [10] it was argued that destructive interference due to T violation reduces Eq. m to the 
approximate bievolution equation Eq. da. The highest degree of accuracy entailed all terms in Eq. m 



0 n 2n 3n 4n Sti 671 


z 

Fig. 6 (Colour online) Plots of \YN-n,n{z)\ = |Tjv-n,Ti(-jy^) | for different values of N verifying that the 
position of the subsidiary maxima are relatively fixed at z = (2m -|- l)7r. 
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vanishing except for n = 0 and n = N which imposes the condition on the total time 


Nt » 


10-13 

VT 


s , 


whereas the less accurate form of the bievolution equation 


(39) 


\mr)) 


^N-n,n + ^ Sn-t 


.n^O 




IV'o) 


is satisfied by the less-stringent condition 

Nt > 10-11' s . 


(40) 


(41) 


However, subtle details of the interference were overlooked. The present work has revealed the in¬ 
terference function Ij^-n,n{z) contains subsidiary maxima that were not previously considered. Their 
presence was shown to lead to a new condition Eq. (1571) that must be met for the bievolution equation 
to be valid. The fact that Eq. dSIl) conflicts with Eq. daH) means that the destructive interference 
cannot be so strong as to eliminate all terms in Eq. ([5]) except for n = 0 and n = iV, as previously 
thought. Instead, the best one can achieve is the more approximate from given by Eq. (HUl) . In that 
case the combination of Eq. (HD) and Eq. (1571) gives the range of values of the total time Nt over which 
the bievolution equation is valid as 

IQ-ii' s< Nt 

We found that the range of allowed Nt values is unlikely to extend to the current age of our universe 
as that would require an unreasonably small proportion of T violating particles. 

Nevertheless, we also found that the new condition condition Eq. dSTJ could be avoided if t, 
rather than being fixed at the Planck time, reduces in proportion to If-^/N + I as the number of 
steps N increases. In that case the upper limit to the total time Eq. dsa no longer applies. However, 
further analysis of this new approach is beyond the scope of the present work and will be presented 
elsewhere |46) . 

In conclusion, although we have shown that the effect of destructive interference for fixed r has 
been overstated in previous work |40j . nevertheless, we have also found that the destructive interference 
can be recovered by modifying the method and allowing r to reduce as the number of steps increases. 
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